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-4— > Abstract 

The influence of nonmagnetic central node defect on dynamical properties of regular square-shaped 5x5 segment of magnetic dot 
array under the thermal activation is investigated via computer simulations. Using stochastic Landau-Lifshitz-Gilbert equation we 
simulate hysteresis and relaxation processes. The remarkable quantitative and qualitative differences between magnetic dot arrays 
with nonmagnetic central node defect and magnetic dot arrays without defects have been found. 
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1. Introduction 

The term Magnetic dot array (MDA) [1,2,3] refers to the 
family of nano-scaled monolayer structures consisting of 
the identical magnetic nanoparticlcs, called dots, which are 
periodically ordered on a non-magnetic substrate. MDA 
concept is compelling, partly because of qualitatively new 
properties, that essentially differ from those of the bulk ma- 
terials. Typical for MDA physics is an intricate collective 
behavior. Its understanding may be valuable for condensed 
matter physics, material science and nanoscience. The spe- 
cial properties of MDA systems follow from the interplay 
between intra-dot and inter-dot interactions as well as from 
the interplay of anomalously large surface compared to the 
bulk magnetic energy contributions [4,5]. The MDA prop- 
erties are already utilized in technological applications con- 
cerning magnetic field sensors [6,7,8] and reading heads of 
magnetic-disk data-storage devices [9] . 

The technology of fabrication of MDAs [10,11,12,13] has 
been perfected to an excellent degree in the past few years. 
In practice, however, an occurrence of technological defects 
and local irregularities has still a great influence on all of 
the magnetic properties [14]. In this paper we study how 
the imperfection in form of the single-dot vacancy affects 
remagnetization of MDA under the assistance of thermal 
activation. 

In the series of papers [15,16,17,18,19] the authors focus 
on remagnetization of small segments of arrays by simulat- 



ing elementary models where dots arc treated as interact- 
ing point dipoles. In Ref. [20] we draw the effect of uniaxial 
anisotropy induced by the eccentrically placed node defect. 
Later [21] we have investigated how the central defect af- 
fects the quasi-static zero-temperature remagnetization of 
MDA. From this study we know that: (i) the square lattice 
seems to be more sensitive to defect occurrence than the 
triangular one; (ii) zero-temperature differences of hystere- 
sis loops are not very pronounced despite of the remarkable 
changes in the local arrangement of dots. In addition, the 
study has opened question of defect influence on the relax- 
ation modes. 

In the present study the effect of nonmagnetic central- 
node defect, in further referred as defect only, is recon- 
sidered. However, the additional realistic factor included 
here is the thermal activation. In such case the statistical 
treatment of results is necessary. Therefore, in order to dis- 
cern defect consequences, two distinct MDA arrangements 
-defect-free (DF) and defect-including (DI) are compared. 

2. Model 

Because the statistical simulations need computational 
effort we preferred use of elementary model where each dot 
is described by the point magnetic dipole. This simplifi- 
cation is justified for monodomain isotropic nearly spheri- 
cal ferromagnetic particles separated by a sufficient lattice 
spacing several times exceeding a dot diameter. Addition- 



ally, we focus on the small MDA samples where dots are 
placed on square LxL lattice. The magnetic state of i-th 
dot is described by the effective rescaled 3d magnetic mo- 
ment m.i normalized as |rrij| = 1. The inter-dot interactions 
are assumed to be dipolar and described by the effective 
field 
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where r,-j is the distance between z-th and j-th dot in 
lattice-spacing units a. The field is measured in the Ho = 
yA/ s (47ra 3 ) -1 units including the dot volume V and satu- 
rated magnetization M s . To study the DI arrangement we 
assume that defect is represented by the zero magnetic mo- 
ment. The dynamics of magnetic moments is described by 
the stochastic Landau-Lifshitz- Gilbert equation [22] 

= -m 4 x hf - ami X (m, x hf) , (2) 

dr 

where a is the dimensionless damping parameter, r is the 
time in to = 47ra 3 [7(l + a 2 )V M^ 1 units and 7 is the gyro- 
magnetic ratio. The effective field hf includes the dipolar 
field hf p , the external field h oxt , and the Langcvin ther- 
mal field hf 1 . Finally hf = hf p + h cxt + hf 1 . The random 
thermal field h* h is defined by averages [23] 

(ft&(r)>=0, (3a) 
(h^hf^r')) = 2DS^ n S(r - t') , (3b) 

where £,77 G {x,y,z} and i, j are the site indexes; D is 
the noise amplitude. According to fluctuation-dissipation 
relation [22,23], the factor D is linked to the temperature T 

- T „ ^ V 2 M 2 
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where To is the characteristic temperature scale. 

3. Hysteresis 

In all of our numerical experiments we set L = 5 for 
MDA. The main empirical argument supporting our choice 
a = 0.1 is that ferrite nanoparticles have the mean value 
of a of such order [24]. Let us first study the remagneti- 
zation in DF and DI systems in the time varying exter- 
nal magnetic field applied in parallel to one of the main 
MDA axes h GXt (r) = (^ xt (r), 0, 0). For simulated geome- 
try similarly as in the hypothetical experimental setup the 
quantity of interest is the magnetization projection M x = 

1 L2 

— mi ■ e x . The component ft,° xt (r) has been cycled 

i=l 

within the bounds — h max < ft- x xt (T) < /i max - The simu- 
lation starts from the nearly saturated state with M x ~ 
1 at /i° xt = /i max , where /i ma x is the bound chosen to 
keep the system of moments saturated. In the rcmagnc- 
tization regime each time-integration step At = 10~ 2 is 
accompanied by the unique change of the external field 
A/i° xt = ±10~ 6 . Both quantities define the sweeping rate 
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v h = (|A/i cxt |/Ar)w = 10" 4 u , where v = H /t . We ob- 
served that numerical results depend on the sweeping rate. 

Numerical integration of Eq. (2) has been performed us- 
ing stochastic predictor-corrector Heun scheme [22]. This 
choice is justified by the fact that in general, the statistical 
error of the scheme can be made arbitrarily small by aver- 
aging over the number of stochastic paths [25]. We tested 
numerical stability of scheme for two different integration 
steps. In order to fix both A/i x and Ar have been 
rescaled by 1 /2. The test confirmed the invariance of statis- 
tical results with respect to the rescaling. Under the condi- 
tions of thermal activation we recorded and treated assem- 
bly of 600 independent loops. The treatment assumes av- 
eraging of the magnetization data conditioned by A/i° xt < 
and A/i x xt > 0, respectively. In Fig.l we show averaged 
loops constructed for T = 0.01T , 0.1T and 0.5T . We see 
that loops differ for DI MDA and DF MDA variants. The 
difference confirms some anomalous impact of the local as- 
symetry of the couplings broken by defect. Clearly, by in- 
creasing the temperature the hysteresis vanishes due to re- 
duction of the impact of irreversible processes due to config- 
urations separated by energy barriers. When assume stan- 
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Fig. 1. The averaged hysteresis loops obtained for DF (left panel) 
and DI (right panel) MDAs at temperatures (a) T = O.OTTo, 
(b) T = 0.1T , (c) T = 0.5T . 

dardly that the amount of free energy responsible for the ir- 
reversible processes is proportional to the area A of hystere- 
sis loop. According Fig. 2 the thermally induced hysteresis 
reveals maximum of the ratio Adi/^4df at T = 0.25Trj. It is 
clear that for sufficiently high temperatures the differences 
between DI and DF systems vanish. Since the remagnetiza- 
tion process of MDAs is very inhomogencous we analyzed 
the local mean magnetic hysteresis for nearest neighbor- 
hood of central node, corners and middle nodes of MDA 
edges. The local loops differ among each other. From their 
comparison it follows that defect affects not only hysteresis 
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Fig. 2. The ratio of the areas of DI and DF hysteresis loops as a 
function of the temperature. 

of its nearest neighbor dots, but also hysteresis at corners 
and edges, that clearly follows from the long-range nature 
of magnetostatic couplings. 

One can expect that averaging eliminates some kind 
of the relevant defect-sensitive information. Therefore, it 
seems valuable to pay attention to the fluctuations around 
the averaged loops. To do this we analyzed the noise formed 
by returns AM X = M x (r + lOOAr) - M x (r) accumulated 
during reversals. Also this data have been treated sepa- 
rately for A/i x xt > and A/i x xt < sweeps. Evidently, the 
probability density functions of returns exhibit defect sen- 
sitivity. They has been quantitatively characterized by the 
parameters of leptocurticity (AM*)/3(AAf^) 2 and skew- 
ness (|AM X | 3 )/(|AM X |) 3 . Their numerical values are listed 
in Table. 1. The values clearly indicate that magnetization 
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0.01 


66.52 


41.31 


2.28 


0.81 


0.10 


4.41 


6.91 


0.05 


0.06 


0.50 


1.93 


1.93 


0.04 


0.03 



Table 1 

Statistical characteristics of the magnetization returns. The leptocur- 
ticity and skewness obtained from fluctuations of hysteresis at dif- 
ferent temperatures. 

noise is strongly non-gaussian. We see that temperature 
affects the characteristics in common manner, they de- 
crease. However, the remarkable differences between DF 
and DI systems exist again. 

The reversal magnetization paths uncovered complex 
collective inner behavior of dot moments. In addition the 
thermal fluctuations cause that single rcmagnctization 
events qualitatively differ between each other. We observed 
that intervals with smoothly varying magnetization, where 
the quasi-coherent rotation of moments prevails, are bro- 
ken by the irreversible jumps. To understand redundant 
but essential attributes of the reversal statistics we ex- 
ploited the abilities of artificial neural networks. The most 
efficient for our purposes seems to be the usage of unsu- 
pervised self-organizing maps (SOM) [26,27] that allow us 
to extract several representative paths from the largely 
redundant assembly In Fig. 3 we show the representative 



paths (T = O.OlTo) extracted by SOM. For certain inner 
parameters of SOM the classification yields two represen- 
tative groups of loops (classification is done for DF and DI 
cases separately). The quantitative output of this analysis 
is that the averaged loops from Fig. 1(a) can be under- 
stood as a mixture of 89% of (al) constituent with small 
addition of 11% of (a2) in DF case, and 65% of (bl) with 
35% contribution of (b2) in DI case. This classification ex- 
plains why the patterns (al), (bl) are much more similar 
to the averaged loops than less probable types (a2), (b2). 
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Fig. 3. The hysteresis loops for T = O.OlTo revealed by SOM network 
for DF and DI array. The boldface branches belong to A/i cxt < 0. 

By analyzing the revealed loops we conclude that two 
massive magnetization jumps contribute to the reversal. 
Deeper insight to the mechanism of their nucleation is of- 
fered by the configuration of moments. The sequence of 
snapshots that correspond to A/i x xt < regime is shown 
in Fig. 4. We see that the defect evidently supports the nu- 
cleation of the intermediate inter-dot leaf. Without defect 
the noncolinear antifcrromagnctic order is prefered. 
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Fig. 4. The configuration of magnetic moments along the hysteresis 
path that belongs to representative loop types obtained for decreasing 
ft x xt equal to (1) 1.5, (2) 0.0, (3) -1.0, (4) -3.0. 



4. Relaxation process 

The additional principle of defect detection is provided 
by relaxation process in zero external field starting from 
M x = 1 state. The relaxation curves for DF and DI sys- 
tems arc plotted in Fig. 5. For T = O.OlTo the qualitative 
differences in the relaxation of magnetization become clear. 
The most anomalous aspect of relaxation of DI is the pro- 
nounced peak. However, DF MDA does not exhibit this 
feature. The configurations show that in DF case the pair 
of antiparallel vortices is formed in contrast to inter-dot 
flower observed in DI case. As for the hysteresis, when the 
energy barriers are over bridged by thermal fluctuations, 
the transition to relaxation mode without peak is expected. 
This scenario is confirmed by the inset of Fig. 5, where the 
nearly exponential relaxation occurs without substantial 
marks of defect sensitivity. 
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Fig. 5. The averaged magnetization relaxation obtained for DF 
and DI systems at temperatures T = 0.01Tb. Inset corresponds to 
T = 0.5Tb . The final curves are averages obtained from 1000 stochas- 
tic relaxation events. 

5. Conclusion and discussion 

The simulation statistical study discusses the theoretical 
grounds for defect detection in MDA. It has been shown, 
that that the presence of non-magnetic central node defect 
in a a small segment of MDA has non-negligible influence 
on dynamical properties of MDA, concretely magnetic hys- 
teresis and relaxation. Further, it shows that the including 
the thermal activation is not only the additional factor of 
realistic simulation, but also the factor that can substan- 
tially enhance the chance to identify defect. 

Finally, we are concluding with the believe that our 
present results will affect the experimental claims in this 
direction. 
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